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QUANTITATIVE ESTIMATES OF THE METACHROMASIA REACTION 
OF VOLUTIN GRANULES OF YEAST USING NEURAL NETWORKS 


Abstract. The metachromatic coloration of volutinous granules of the yeast Saccharomyces cerevisiae is one of 
the indicators of the influence of sharp geomagnetic field (GMF) perturbations. The metachromasia reaction is based on 
the aggregation of dye molecules in interaction with inorganic polyphosphates, which are components of volutinous 
granules. To determine the characteristics of the geomagnetic field that cause the appearance of different colors of the 
metachromasia reaction, it is necessary to simultaneously monitor this reaction and changes in the GMF. High-quality 
monitoring is possible with rapid automated counting of cells with all possible color changes during the metachromasia 
reaction. The aim of the work was to develop a neural network architecture for recognizing and quantifying color changes 
and heterogeneity in real time during monitoring of the metachromasia reaction of volutinous granules of the yeast S. 
cerevisiae, which is necessary for further determining their correlations with changes in the geomagnetic field of different 
intensities. A program based on a nonrecursive labeling algorithm was created to count the number of cells in the study 
groups. In the course of the work, the software of two neural network architectures was compared to determine the best 
results in recognizing and quantifying yeast cells with different colors during the volutinous granule metachromasia 
reaction. It was determined that the Unet architecture type coped with the tasks of cell classification and segmentation 
much more efficiently than the Inception v3 architecture. The average relative error for automatic recognition of all cell 
groups was 3.85%, and the maximum relative error was 4.56%. The performance of the neural network was 89.9% when 
detecting cell segmentation and 86.4% when detecting color differences in the metachromasia reaction. 

Keywords: metachromasy, geomagnetic field, volutin granules, neural networks, yeast. 


Introduction concentration and chromotropic effects that 

Due to significant climate change and the depend on the electromagnetic background 
impact of space weather on biota, it is an urgent [11, 12]. 
issue to identify biological systems that are An important task related to the study of 
sensitive to changes in the environment [1-5]. the number of cells with changes in the color 
This can be useful in the development of of the metachromatic reaction of volutinous 
biosensor devices that are sensitive to sudden granules during geomagnetic field 
changes in geophysical factors [6, 7]. It has perturbations of different strengths is to 
been established that the color change of the develop methods for rapidly accounting for 
metachromatic reaction of volutinous granules such changes. The manifestation of the color of 
of the yeast Saccharomyces cerevisiae is a the metachromatic reaction cannot be 
potential bioindicator of geomagnetic field measured by the colorimetric method, since 
perturbations [4, 8-10]. In this case, the aqueous dye solutions do not obey Beer's law 
phenomenon of metachromasia is manifested [5]. Technologies for manual visualization of 
in a characteristic color change (different from changes in cells are still widely used in cell 
the color of the "methylene blue" dye) in the counting methods. However, they have a 
process of staining yeast volutinous granules. number of disadvantages, such as low 
It has been shown that the main colors of sensitivity and specificity, low reproducibility 
metachromasia are largely due to the of results, and a long time to detect and 
aggregation of dyes, which is induced by evaluate changes in cells [5, 13]. Such methods 
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of assessing changes in the color of the 
metachromatic reaction are not capable of fast 
and automatic data processing during 
monitoring studies, which are necessary to 
determine the correlations between changes in 
the color of the metachromatic reaction and 
geomagnetic disturbances [4, 5, 13, 14]. For 
high-quality monitoring, it is necessary to 
quickly automate the counting of cells with all 
possible color changes during the 
metachromatic reaction. 

It is known that the reorganization of 
electronic interactions within the aggregates 
leads to a shift in the resonant frequency, 
which is manifested in 3 types of color changes 
in the metachromatic reaction [15].When 
staining, we distinguish between:- a- 
metachromasia (non-metachromatic color), B- 
metachromasia (weak purple color, 
intermediate between blue and red), and y- 
metachromasia (strong positive red color) [4, 
15].These colors indicate whether the dye is in 
one of the following forms: a dye monomer (a- 
metachromatics), a short polymeric dye of two 
and three or more components (f- 
metachromatics), or a_ structurally more 
complex polymeric dye (y-metachromatics) [5, 
12, 17].However, yeast volutin grains often 
contain segments of heterogeneity where the 
color of the metachromasia reaction changes 
[18, 19], which is manifested by the 
appearance of red areas within still blue volutin 
granules [5, 12]. In addition, there are a certain 
number of cells in which no color change was 
observed during the metachromasia reaction 
due to the absence of a reaction, and they must 
also be counted. Therefore, for qualitative 
monitoring studies, it is necessary to recognize 
and quantify 5 types of changes in cells during 
the metachromasia reaction: 3 types of cell 
color changes, cells with heterogeneity, and 
cells with no reaction [20, 21]. 

The use of neural networks with certain 
training makes it possible to automatically and 
in real time assess changes in color in cells and 
take them into account during _ the 
metachromasia reaction [22-25]. Such an 
approach to assessing changes in the number 
of cells with different colors is an effective way 
to determine the degree of influence of 
geomagnetic field perturbations of different 
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intensities on the metachromatic reaction of 
volutin granules of the yeast S. cerevisiae. 

The aim of this work was to develop a 
neural network architecture for recognizing 
and quantifying color changes and 
heterogeneity in real time when monitoring the 
metachromasia reaction of S. cerevisiae 
volutinous yeast granules of UCM Y-517 and 
determining their correlation with changes in 
the geomagnetic field of different intensities. 


Materials and methods 

The dataset of twenty photographs of 
cells S. cerevisiae UCM Y-517 colored in 
methylene blue (by Lefler) [5] is used in the 
paper. A photo editor, "Corel PHOTOPAINT 
X6," was utilized to process photos. Taking 
into account the peculiarities of the dataset 
obtained, the methods of pre-processing 
images for segmentation of cells with volutin 
granules according to different types of color 
metachromasia response and cell classification 
according to their classes of presence and 
absence of metachromasia, respectively, have 
to be determined [26, 27]. In order to count the 
number of cells with a particular characteristic 
(index k in 5 degrees) of the volutin granule 
metachromasia reaction (S), the following 


formula was used: S = nee , where 


D, is the presence of a cell metachromasia 
reaction of the k-type, and Nx — is the number 
of determined cells of a given type. 

When building the NNs, it was necessary 
to choose the network architecture type, select 
learning scales, set the number of inputs and 
transfer functions of the used neurons and their 
interconnection, and, finally, determine the 
outputs of the particular network [28 - 30]. For 
this purpose, a convolutional NN was selected 
using the Python programming language and 
the Keras open NN library [31 - 33]. The 
Inception v3 and Unet models were applied 
when modeling the architecture of the 
particular NN. Those were earlier developed 
for the problems of classification and 
segmentation of objects in the image [34,35] 
and for cell segmentation in medical images 
[36, 37], respectively. 

To solve the tasks of cell segmentation 
and their classification into two classes of 
presence and absence of metachromasia, the 
program codes for the NN were proposed 


ISSN 2710 - 1673 Artificial Intelligence 2024 Ne 2 


according to the two NN architectures. Then 
the NN was trained using the obtained data. In 
addition, the effectiveness of trained NNs was 
tested using test data. Then, as a first step, cells 
were classified and segmented by applying a 
trained NN. The number of cells was 
calculated using a non-recursive markup 
algorithm as a second step [38, 39]. The 
statistical calculations of the color 
characteristics of yeast cells with a detected 
metachromatic reaction and cells without this 
reaction or of the NN efficiency for two NN 
architecture types were performed separately 
using the MS Excel 2010 program. 

Further, it was necessary to build a deep 
convolutional NN [40 - 42] designed 
specifically for processing input images of 
volutin yeast granules. The main idea of this 
process was to supply the usual network with 
additional sequential layers where increment 
sampling operators replace union operators. 
Hence, these sequential convolutional layers 
increase the output resolution due to increased 
digitization [43, 44]. The high-resolution 
functions were joined with digitized outputs 
via feedback to perform cell localization. It 
enables training a sequential convolutional 
layer to get a more accurate result using 
additional information [45, 46]. 

The Inception architecture was used to 
perform the mentioned operation. It consists of 
two main blocks. The first block determines 
the NN features and plays the role of an 
extractor of functions. The block ensures 
compliance with the templates by applying 
convolution filtering operations to fulfill this 
function. The first aggregate layer filters the 
image using multiple kerners and returns 
"object maps, " normalized with the activation 
function or resized. Recognition of objects on 
the input layer was represented as a three- 
dimensional grid, the dimensions of which 
depend on the input image, according to the 
formula: 

Ih=W-H-D, 


where J, is the input NN layer size, W and H 
are the widths and the height of the input 
image, respectively, D is the number of color 
channels of the image [47]. A maximizing 
aggregation layer was selected with 22 filters 
used in increments of 2, which reduces the 
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sampling of each section of the input depth 
twice for both width and height, discarding 
75% of the excitations. It is worth stressing that 
each operation takes the top acts over four 
numbers, keeping the depth size unchanged 
[48]. 

The NN architecture of GoogleNet 
Inception v3, trained on about a million images 
and about a thousand object classes, was taught 
on our dataset using training with a teacher [37, 
42, 49, 50]. The transfer learning method is 
used for the NN. This method allows for 
retraining the last level of the existing model in 
order to significantly reduce the learning time 
and size of the necessary data chosen. This 
paper uses the Inception v3 model trained for 
the ImageNet Large Visual Recognition 
Challenge on 2012 [51] data as a pre-trained 
model. 

The NN based on the Unet architecture 
consists of a narrow path corresponding to the 
convolutional NN's type and architecture. That 
is represented by reapplying two 3x3 
convolutions, followed by a layer of offcut line 
nodes (ReLU — Rectified Linear Unit) and the 
2x2 aggregation operation in increments of 2 
to reduce the resolution [42]. ReLU uses the 
nonviolent transmitting function f(x) = max (0, 
x) for these operations to fulfill this step. The 
ReLU_ correctional layer enhances the 
nonlinear properties of the decision-making 
function and NN as a whole without affecting 
the recipe fields of the convolutional layer. We 
have changed this architecture to work with a 
few learning images and give more accurate 
segmentations. 

Because of very little training data, we 
have used artificial data magnification, 
applying deformations to existing learning 
images. It allows the NN to study the 
invariance of such deformations, which is 
essential since, in biomedical segmentation, 
deformations are the most common variations 
in tissues and cells. The importance of 
increasing the data for invariance teaching in 
uncontrolled learning of properties was shown 
in [42]. 


Results 

In this paper, the analysis of the dataset 
of 20 images was carried out. The dataset 
features were established, considering the 
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proper choice of methods for pre-processing of 
images. A training sample of photographs was 
prepared after pre-processing the initial dataset 
of images. To reduce the requirements for 
computational power, all images were 
transformed from a size of 1600x1200 pixels 
to that 224x224 pixels. For the NN architecture 
types based on the Inception v3 and Unet 
models training masks of individual segmented 
cells were created using the Corel 
PHOTOPAINT X6 photo editor (Fig. 1). The 


Then, the image portions of a lower than 
average brightness were determined, those 
being surrounded by pixels of a higher than 
average brightness. _—_ Binarization _for 
segmentation of each cell was applied (Fig. 
2b). Further to that; the NN is going to carry 
out the classification of segmented cells 
denoting the corresponding class in color: 
"presence of metachromasia" - red, "absence of 
metachromasia" — blue (Fig. 2c). 


NN training for cell segmentation and ins wns st AHR ets Pr) We 
classification according to two features of the a AS Nie’ eign its Ba 
cell reaction, namely presence and absence of THe HANG ANes' Weis tyes 
: ; ayiN ° Si v5 o2 jel Ne. { hy} 
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The analysis of the metachromasia is Wd inert IE N ONG, aie 
reaction of yeast cells S. cerevisiae. Python Wy) Wee vo Wels yon SAN 
: V2 <NG a TN 7A en 
programming language created the NN based ely bir} Nis “64 seniliaises 
on the Inception v3 architecture. The NN code mUhes aN ihr Ht a Hye! 
was developed to recognize cell segmentation “as SENS Nahe Nisan 
according to metachromasia's presence or 4 aH i aS y Us 
absence. Then the NN training with testing the ah HH Z ~~ ais ay 155: 
effectiveness of the trained NN was carried out \ gil wal HT M4 Messy: we asst 


using test data. An average brightness for the 
image portion was determined first (Fig. 2a). 


Figure |. Training mask for teaching the NN 
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Figure 2. NN (based on Inception v3 architecture) input and output images: a) input image, b) image with segmented 


cells, c) image with classified cells 


Figure 3. The image with identified cells 
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After the trained NN already segmented 
and classified cells, the number of cells with a 
metachromasia reaction can be calculated 
using a non-recursive markup algorithm. The 
NN based on the Inception v3 architecture 
analyzes the image and produces results of the 
analysis of the neuron in the form of red 
rectangles by which the detected cells are 
visualized (Fig. 3). On top of that, an analysis 
of the NN effectiveness was carried out using 
methods of mathematical statistics. As visible 
in Fig.3, the NN based on the Inception v3 
architecture could not isolate all cells with 
existing features. The program effectively 
performs cell segmentation but makes errors at 
the next stage of cell classification. When the 
volutin granules of cells do not have a clear 
blue color, the NN classifies them incorrectly 
as "metachromasia." This type of architecture 
needs sufficient detection accuracy. It requires 
significant computing resources (highly time- 
consuming) since many neural layers are used. 
The NN efficiency was found to be 80.3%, 
which means an insufficiently accurate cell 
separation and the necessity to increase 
training and work efficiency data to receive 
satisfactory results [52, 53]. Naturally, to 
continue NN training and improve the 
methodology of pre-processing of images, 
particularly reducing noise, is the way to 
distinguish an image according to the actual 
size of pixels instead of approximating it as a 
whole. 

To address these difficulties, it was 
decided to conduct the NN (based on the Unet 
architecture) training using training images 
and training masks of segmented cells. The 
main difference between the Unet architecture 
and the Inception v3 architecture is in the way 
of teaching and visualization result. The Unet 
NN can carry out adequate cell segmentation 
in the images and highlight the silhouettes of 
objects by binarization. The data obtained were 
presented in the dataset of 20 images. Initially, 
a dataset with a volutin — granules 
metachromasia reaction in yeast cells S. 
cerevisiae was separated into two sets - the 
training set and the set for verification. 
Preliminary processing and expansion of the 
image, which included standardization and 
scaling of the input image to improve the 
dataset transfer quality, was carried out. 
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Training masks of segmented cells were 
created using the Corel PHOTOPAINT X6 
photo editor. 

In preliminary data preparation, the 
images' main shortcomings were leveled. 
Those are the blurred contours and boundaries 
of cells, the presence of dye spots in the 
background, the lack of precise color gradation 
depending on the type of metachromasia 
reaction, differences between cell sizes, and 
too large cell size. These shortcomings 
increase noise and reduce the data processing 
efficiency. In order to expand the dataset, the 
images were adapted to the size of the NN 
model; mask fragments of 224x224 pixels 
were created from a size of 1200x1600 pixels. 
To improve the model's accuracy, images of 
individual cells were added to the dataset to 
move from the overall rating format to the 
specific one. 

As noted above, two datasets were used 
in the paper: a training set and a set for 
verification. A dataset fulfilled the verification 
set's role of adjusting the model's 
hyperparameters for training. 

When creating the NN based on Unet, 
the Python 3.7 language was used, and the NN 
code was developed. It improved the NN 
efficiency. The code distinguished three stages 
in the process of analyzing the metachromasia 
reaction. At first, the segmentation of yeast 
cells in the image was carried out by the NN, 
followed by the classification of cells by the 
degree of the color of the reaction of the 
volutin granules of each cell in the image. The 
degree of metachromasia reaction for the entire 
image was then estimated along with checking 
the total number of cells and the number of 
cells for which the metachromasia reaction 
was present. In other words, to classify cells 
the image with segmented cells without 
classification was obtained before obtaining 
the image with cells classified as "presence of 
metachromasia" and finally that of "absence of 
metachromasia". 

Counting the number of cells. To 
calculate the number of cells in 5 considered 
groups, the program for counting the number 
of cells has been created. It is based on a non- 
recursive markup algorithm. This algorithm 
works via the utilization of a corner-like shape 
so-called the ABC mask [34, 54] (Fig. 4). To 
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find out specific color components, the passage 
of this mask through the images of cells takes 
place from left to right and from top to bottom. 


Figure 4. ABC mask and direction of sequential image 
(its green part) scanning 


Five possible positions of the mask from 
top to bottom on color images of 
metachromasia analyzed are shown in Fig. 4: 

e position 0, when all three components of 
the mask are not marked (color reaction of 
metachromasia 1s absent) — in this case, a pixel 
is Skipped (no mark is created); 

e position 1, when only the color 
component A (a-metachromasia) is marked — 
in this case, a mark | is created; 

e position 2, when only the component B 
(B-metachromasia) is marked — in this case, 
pixel A is marked as 2, which is located in B; 

e position 3, when only the component C 
(y-metachromasia) is marked — in this case, 
pixel A is marked as 3, which is located in C; 

e position 4 means that the labels (object 
numbers) B and C are equivalent, and pixel A 
can be labeled as B or C. If the value B is not 
equal to the value, C, all already processed 
pixels labeled as C are relabeled in B and 
renumbered. 

Then the number of cells, namely the 
total and number of cells for each of the three 
types of metachromasia reaction, is counted 
using a non-recursive markup algorithm. 

After the dataset analysis, the network 
segments individual yeast cells, defines the 
area of interest where the metachromasia 
reaction can be detected, and counts cells with 
an apparent metachromasia reaction. The 
method of automatic semantic segmentation 
was used for image segmentation. Using this 
approach, the model works with a set of objects 
in which each pixel is labeled. That is used in 
further analysis to be carried out (Fig. 5). 
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Figure 5. Comparison of the input test image (left) with 
the output image given by the NN based on the Unet 
architecture 


A comparison of the cell color 
background when each pixel is assigned a 
number was used in this analysis. These 
numbers helped indicate the area in which the 
phenomenon of metachromasia was detected. 
Such zones stood out with yellow marks, 
significantly narrowing the range of reactions 
to be searched for the color. The white marks 
determined the most contrasting pixel inside 
the range of yellow marks, which is necessary 
for further counting the number of cells with 
metachromasia. Attaching marks in the form of 
dots or crosses avoided blurring the data with 
a rectangular frame superimposed on each 
other, thereby reducing the result's accuracy 
(see Fig. 3). 

The Inception v3 NN model was used to 
count the total number of cells and also to 
count separately both the number of cells for 
the present (3 types of images corresponding to 
1 — a-metachromasia; 2 - B-metachromasia; 3 - 
y -metachromasia, as mentioned above) and 
absent metachromasia reaction of volutin 
granules of yeast S. cerevisiae. Results for 
manual and automatic methods are presented 
in Table 1. 
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Table 1. Manual and automatic counting of the number of cells in datasets with (three types of metachromasia reaction) 
and without metachromasia 


Type of Total number of counted Number of cells with Number of cells with 
metechro- cells "metachromasia" "no metachromasia”" 
pera Manual Automatic Manual Automatic Manual Automatic 
1 363411 374+10 16 20 347 354 
2 474422 496+23 26 28 448 473 
3 324413 337414 34 4] 290 296 


The average relative error value for 
automatic recognition of these three groups of 
cells (manifesting three types of 
metachromasia reaction) is 3.85%, the 
maximum relative error being 4.56%. The 
Inception v3 NN model has coped with the 
recognition of cell contours. However, 
difficulties have arisen with recognizing cells 
manifesting a metachromasia response. It was 


influenced by the insufficient image quality, 
insufficient quality of their processing, and the 
amount of data in datasets. The training of the 
Unet NN model made it possible to determine 
the number of segments with small isles of the 
metachromasia reaction and the number of 
segmented cells without such a reaction in 
segments (Table 2). 


Table 2. The Unet NN results of determination of segments with three types of metachromasia reaction in cells 


Number of segments with ‘ 
3 : Number of segments without 
Type of metachromasia (3 types of reaction) rae 
bat : metachromasia in cells 
metachromasia image in cells 
ee Manual | Automatic Ferme Manual | Automatic Relea 
error, % error, % 
1 0 0 - 27 30 10 
5 9 20 16 18 11 
3 36 39 7,6 43 47 9,3 
Discussion groups of cells was revealed to be 3.85%, and 


As can be seen from Table 1, in cells 
with type 1 metachromasia reaction, all 
detected segments were free of 
metachromasia. In the case of the second type 
of reaction, a few segments were found with 
the metachromasia reaction and only a few 
segments without metachromasia in the images 
of segmented cells. The most significant 
number of with and without apparent 
metachromasia responses were easily 
identified on images with the third type of 
reaction (Table 1). At the same time, some of 
the identified segments did not have 
information about the color change of the 
metachromasia reaction. However, it should be 
noted that the number of cells with accurate 
results obtained by registration using a neural 
network is higher than that obtained by manual 
registration [12, 30, 52). 

The study revealed the advantages of the 
automatic method for processing 
metachromasia reaction data. The average 
relative error for automatic recognition of all 
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the maximum relative error was 4.56%. The 
neural network can be used to automate the 
estimation of the number of cells and the 
intensity of their coloration in monitoring 
studies of the impact of geophysical 
disturbances on the metachromatic reaction of 
volutinous yeast granules [4, 5]. 

Training of the neural network model 
allowed us to determine the number of 
segmented cells with the presence of the 
"metachromasia reaction" in the segments and 
the number of segmented cells with the 
absence of such a reaction in the segments 
(Table 2). The average relative error was 
10.1% in the identification of segments in cells 
without a metachromasia reaction and 13.8% 
in the identification of segments with a 
metachromasia reaction (Table 2). The 
developed program using the Inception v3 and 
Unet architectures for detecting and processing 
color changes in the metachromatic reaction 
showed that the model of such a convolutional 
neural network effectively copes with the tasks 
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of cell segmentation recognition at the level of 
89.9%, and when detecting the types of 
metachromatic reaction in the segments, the 
recognition rate was 86.2%. This indicates the 
effectiveness of the NN at a sufficient level of 
>80% [55] to detect color changes in the 
metachromatic reaction using the applied 
specific neural network training. 

Moreover, the great advantage of the 
automatic method of counting the cells is its 
higher accuracy and rate of calculations, which 
is extremely necessary for monitoring studies 
of cell transformations and in_ further 
application of data obtained for 
interdisciplinary research [5, 19, 22]. 

The neural network coped with the 
recognition of cell contours, but there were 
difficulties with the recognition of cells with 
the absence and presence of a metachromasia 
reaction. This was due to poor image quality, 
insufficient processing, and the amount of data 
in the datasets [24, 42, 53]. In the future, the 
training of the neural network will continue 
and the methodology for image quality pre- 
processing will be improved: noise reduction, 
which will significantly increase the efficiency 
of the neural network, as well as the need to cut 
the image by a certain pixel size, rather than 
bringing it to this value as a whole [36, 47, 50]. 


Conclusions 

To solve the problem of classification of 
color changes in volutinous granules and 
segmentation of S. cerevisiae yeast cells during 
the metachromatic reaction, two Python 
programs were developed using _ the 
architecture of artificial neural networks 
Inception v3 and Unet, and it was found that 
the architecture of Inception v3 networks is 
less effective in performing the task of cell 
segmentation, and the neural network based on 
the Unet architecture is optimal. 

The cells classified by color and the 
presence of segmentation were counted using 
a nonrecursive markup algorithm and it was 
found that the relative error of the created 
software does not exceed 4.53% and averages 
3.85%. The use of this neural network 
algorithm allows us to more effectively 
determine the dynamics of the metachromatic 
reaction under the influence of geophysical 
factors in an automatic mode. 
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An accurate and quick assessment of 
both the number of cells and differences in the 
metachromasia reaction color of volutin 
granules of Saccharomyces cerevisiae yeast as 
an effective processing code for 
metachromatic staining is essential not only for 
searching biological problems but also for 
interdisciplinary research. 
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